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A Class of Asymptotic Orbits in the Problem of Three 

Bodies. 

By L. A. H. Waeren. 
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I. Introduction. 

If we have a system consisting of two arbitrary finite spherical bodies 
revolving in circles about their common center of gravity, Lagrange* has shown 
that there are three points on the straight line passing through the centers of 
the finite bodies such that if an infinitesimal body be placed at one of them and 
be projected so as to be instantaneously fixed relatively to the revolving 
system it will always remain fixed relatively to the system. These three 
equilibrium points, as they are called, are separated by the finite bodies, whose 
masses are denoted by (i and l—[i, {0<[i<i). The point beyond the mass fi 
is called (a), that between the finite bodies is called (&), and that beyond the 
mass 1 — (I is called (c). 

In the present paper, it is shown that in the neighborhood of each of the 
equilibrium points (a), (h), (c) there exists a class of orbits, in which the 
infinitesimal body approaches the equilibrium points as t becomes infinite- 
Such orbits are called asymptotic orbits.f 

* Lagrange, " Collected Works," Vol. VI, pp. 229-324 ; Tisserand, " Mecanique Cdleste," Vol. I, Chapter 
VIII; Moulton, "Introduction to Celestial Mechanics" (New Edition), Chapter VIII. 
■fPoincare, "Les M6thodes Nouvelles de la Mecanique Celeste." 
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In § II the equations of motion of the infinitesimal body are given with 
reference to the system of rotating axes, passing through the center of gravity 
of the finite masses, the x-axis coinciding with the line joining their centers. 

In § III there is a detailed explanation of what is meant by asymptotic 
solutions of systems of dilTerential equations, and asymptotic orbits of a 
moving body; and it is shown that in the three-body problem under considera- 
tion the only existing orbits which are asymptotic to any of the equilibrium 
points (a), (b), or (c) lie wholly in the plane of revolution of the finite masses. 

In § IV the equations of the asymptotic orbits are developed as power 
series in exponential functions of the time, and in the following article it is 
shown that these power series expansions are convergent for t sufficiently 
large. 

In § VI an alternative method is given for building up the equations of 
these orbits. 

In § VII there is a discussion of some of the principal properties of these 
asymptotic orbits, their position relative to the rotating axes, and the change 
in their direction of approach to the equilibrium points as fi varies from zero 
toi 

In § VIII it is shown how the orbits can be continued by the method of 
mechanical quadratures beyond the range of convergence of the solutions of 
the differential equations. A special case //=0.02 is discussed in detail, and 
one of the corresponding asymptotic orbits for point (a) is continued by this 
process. 

II. The Differential Equations of Motion of the Infinitesimal Body. 

In the following discussion we consider a system consisting of two finite 
bodies revolving in circles about their common center of mass, and of an in- 
finitesimal body subject to their attractions. Let the constant distance between 
the finite bodies be unity. Denote the masses of the finite bodies by fx and 
1—fi, where 0<|U<|, so that the sum of the masses shall be unity. Choose 
the unit of time so that the gravitational constant F shall be unity. With the 
units so chosen the time of revolution of the finite bodies will be unity. 

Take the origin of coordinates at the center of mass of the finite bodies, 
and refer the motion of the bodies to a system of axes rotating in the plane of 
motion of the finite body, in such a way that the ^-axis always passes through 
the centers of the finite bodies. If ^, rj, ^ denote the coordinates of the in- 
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finitesimal body, then the differential equations of motion for the infinitesimal 
body are :* 



-2-^=^-(i-i")^^^::j^-i" 



dt' dt ~^ ' ^' rl ^ rl 



(1) 



where i^i, 0, 0), {^2, 0, 0) are the coordinates of the bodies 1 — ^u and /w respec- 
tively; ,-i = V'(|=|J^+?+f and r.^VIF^^r+^^Kl 

Let (^0, 0> 0) denote the coordinates of one of the equilibrium points on 
the ^-axis. If, then, by the transformation ^=^0+*') V=y, ^—^, ^6 move the 
origin to one of these points, the equations of motion take the formf 

x"—2y'= {l + 2A)x+ '^B[ — 2x'+y^+s''\ +2C[2af—3xy^—Sxz^] .... 
y" + 2x'={l-A)y + 3Bxy+^Gy[-^a^+y'+z'] :]-•..., 






where 






_ + 1— /i + /i 



5 = 



rf * — rf)* ' 






(2) 



where in the expression for B the upper, middle, or lower signs are to be taken 
according as the orbits in the vicinity of the point (a), (&), or (c) are being 
treated. 

In what follows we shall have to deal chiefly with the first two equations 
of (2), and it will be more convenient to have them in a normal form. 
The linear terms of the first two equations of (2) are 

x"-2y'-{l + 2A)x=Q, I 

y"^2x'-{l- J)^=0. J 

*Moulton, "Introduction to Celestial Mechanics" (New Edition), § 152. 
tMoulton, " Periodic Orbits," pp. 156, 168. 
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The general solution of these equations is 

x=Kie'"^-''+K^e-''^-'*+Kse'"+K^e~''\ 
y^nV^-i[Kie'"^-''—K^e-'''-^']+m[Kse'''—K,e-'"]. 



(4) 



Where ctV— 1, — aV — 1, p, — p are the roots of the bi-quadratic equation 

;i*+ (2-^)^.^(1+2^) (1-^)^0. (5) 

If then we make the transformation 

X =% + M2 + % + 'M4, 



(6) 



(7) 



a;'=<rV— 1(mi— M2)+p(M3— %), 
y =n\/—l{uj,—u<i)+m{us—Ui), 

the first two equations of (3) assume the normal form,* 

2(wcr — np)V — IL^ -, 

+ 2C{2a^—3xy^—3xs^) + ....{ 

- 0/ \ A 3Bxy+lcy{-Ax'+y' + ^) + ... 
2(wp+w(T) 12 

'*^"~*'^^~^'*^~2(m(T-np)V-iL J~ 2{mp+m) [ J ' 

'''^ P"^~ 2{m^-np) L J+ 2{mp+na) I J ' 
<^-p^*+ _^__^ |^"J+ _______ |"| . 

Equation (7) may be put in the form 

< = -pM4+Hf>+flf +HP+ . . . ., 
where Hj^^ denotes all the terms on the right-hand side of the k-ila. equation 
which are of degree r in x, y, z, and therefore of degree r in Wi, u^, u^, Ui. 

III. Asymptotic Orbits Defined and Proof that They all Lie in the xy-Plane. 

It has been shown by Poincaref and Picard| that certain systems of 

differential equations of the form -j^ =Xi{x, t) {i=l, . . . .,n) admit of solu- 

* Moulton, " Periodic Orbits," pp. 161-162. 

f Poincare, "Les Metliodes Nouvelles de la Mecanique Celeste," Vol. I, Chap. VIII. 

t Picard, " Trait6 D'Analyse," Vol. Ill, Chap. VIII, § V. 



(8) 
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tions as power series in Aie^** (^=1> . . ■ ., n), where the Ai are arbitrary con- 
stants, and where the fixed constants \, which are the roots of the characteristic 
equation, are called the "characteristic exponents." It has been shown 
further* that if there are h of the \ (*=1, ....,«) which are represented by 
k points on the complex plane all of which lie on the same side of a straight 
line passing through the origin, and which are such that none of the relations 

h 

^Pj^j — ^i=0 (^=1, . . . ., n) holds for any positive integral values of the Pj 

^=' k 

such that 2 Pj^2, then the solutions as power series in A^e^'* {i=l, . . . ., k) 

will be convergent for | Aie'"' | sufficiently small. In particular, if we put equal 
to zero the Ai corresponding to those ?^i whose real parts are zero or positive, 
the solutions as power series in the remaining A^e^'* will be convergent for all 
values of t which are sufficiently great ; and if the Ai involved in these latter 
expansions are taken sufficiently small the convergence will hold for all values 
of t from ^ = to t = oo . 

Again, if we build up solutions as power series in those A^e^'' where 
the real parts of ?,j are positive, the exponentials e'"'' approach zero as t ap- 
proaches — CO . Such expansions will be convergent for i sufficiently large and 
negative; and if the Aj are sufficiently small they will be convergent for all 
negative values of t. Such solutions are said to be " asymptotic " to the solu- 
tions obtained by putting all the Ai (i=l, . . . ., n) equal to zero, and they are 
called " Asymptotic Solutions of the System of Differential Equations." 

In the problem under consideration we shall show that the differential 
equations of motion of the infinitesimal body have asymptotic solutions such as 
have just been described. We shall see also that the infinitesimal body, mov- 
ing in an orbit defined by one of these solutions, will approach asymptotically 
one of the equilibrium points (a), (6), or (c) as t becomes infinitely great. 

We proceed to show, first of all, that all orbits which are asymptotic to 
one of the equilibrium points (a), (&), or (c) lie entirely in the plane of revo- 
lution of the finite bodies, thdt is, in the a;«/-plane. 

If in the equations of motion (2), we make the transformation x=xs, 
y=ye, s=^ss, we obtain, on dividing through by e, 



x"-2y'-{l + 2A)x = eX,{x\y',z')+B'X,ia>',xf,xs')+s'X^i. . . .) + 
y"+2x'-(l- A)x=6yY,{x) +e^yYs{x% y\ s') +e'yl\{. ...) + .... 
s"+ Az = ezZ^ix) + e'zZs {x^ y\ s^) -\- izZ, (....) + .... 



(9) 



* Poincar6, "Lea Methodes Nouvellea de la Micanique Celeste," Vol. I, Chap. VIII, § 105. 
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Since the right-hand members of (2) converge, so also will the right members 
of (9) converge for all s (0<f <1). It follows, therefore, that the equations 
(9) have a solution of the form 



x= -^xAt)^, 

,=0 

00 

z = 2 zAt)eK 



(10) 



In order that z may be zero for all t sufficiently large, and for all values of e, 
it follows that Zj{t) (i=0, . . . . oo ) must each be zero for all t sufficiently large. 
On substituting (10) in (9) and equating coefficients of like powers of e on 
both sides of the resulting equations, we obtain sets of differential equations 
from which the values of Xj, «/,-, z^ (i=0, . . . .oo ) can be obtained sequentially. 
From the terms independent of e, we have 

<-22/;- (1 + 2^)0=0 = 0, 

y'^ + 2x,-{l- A)y,^Q, I (11) 

4' + Az, =0. 

The general solution of equations (11) is 

^,=nV^[^f>e^^-"— ^fe-^^-^J+wCATe"'— ^fe-"'], ■ (12) 
2o = cf ) cos VI t + cf ■ sin VZ t. 

In order that Za(t) shall approach zero for t infinite, we see that cf^ = cf^ — 0, 
and, therefore 

«o(0=0. (13) 

When we equate the terms in the first power of e, we get 

x'^-'iy':-{l + 2A)x^=X,{a?„ yl, zl), 

y[' + 2x['-{l- A)y^=yoY,{Xo), (14) 

z'l + Asi=SoZ2{xo). 

If we substitute the value Zo=0 in (14) the third equation becomes z'i+Azi=0, 
which is of the same form as the third equation in (11). In order that Zi{t) 
shall approach zero for t infinite we see, therefore, that Zi{t) ^0. Similarly, we 
can show that Z2{t)^0. Suppose we have proved sequentially by this method 
that s^{t)^0, for j=0, 1, . . . ., n. Since the right-hand member of the third 
equation of (9) carries as a factor, and since the factor s has been removed 
from this equation, it follows that the right member of the differential equation 
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which defines s„+i will consist of terms which carry as a factor one or more of 
the Sj{t) (i=0, 1, . . . ., n) but will not contain any Zj, j>n. Hence, the right 
member of this equation is zero, and, therefore, s„+i(0 ^^0- We see, therefore, 
that 0,(0 =0 for i = 0, 1, . . . . oo , from which it follows that z{t) =0. In order, 
therefore, that the infinitesimal body shall come to rest at one of the equilibrium 
points (a), (&), (c) its whole orbit must be in the xy-plane. 

IV. Formal Construction of the Solutions. 

We proceed to show that equations (7) admit of solutions defining asymp- 
totic orbits in the a;i/-plane. These will be found as power series in e""*, and 
6+"*, the former convergent for t sufficiently large and positive, the latter for t 
sufficiently large and negative. 

(A) Solutions as Power Series in e""*. 

If, in equations (8), we make a transformation on the independent variable 
by putting w = e""* these equations take the form 

It is required to find solutions of (15) as power series in w, convergent for w 
sufficiently small. 

By Maclaurin's expansion 

u,{c.)=^u,{0)+c.(^) +4r(^) +••••' (* = 1, •••m4). (16) 

Since the body is to be at rest at the origin at 0=0, it follows from (6) that 
M.(0)=0, {i = l, . . . ., 4). By repeated differentiation of equations (15) with 
regard to a, and putting o = 0, we can build up the coefficients of the successive 
powers of o in the expansion (16). Since 

X = lll + U2 + Us + Ui, I 

y=n^/'IIl(^Ui—U2)+m{us—Ui), J 



(15) 
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it follows that the coefficients of successive powers of o in the expansions for 
X and y will thus be known also. 

The first equation of (15) written at length is 

9% / T r 3mB / o 2 r 2\ '^-^ 



oui , — - r 3mB / o 2 , 2\ ^^ / x1 

— pa-^ = 0V_li(.+ -==. (— 2a;^+i/^)— -^, ; r- ixn) 

+ r ^^ ^ (2a;«-3V) - ,, ^^' , {-4:x'y+f)'] 

+ (terms of higher degree in x and y). 
On differentiating this with regard to a, we get 

— p ^ — ptt ^^ = (T V — 1 ^- 

dw dw dtt 

r 3OT.g / ^^dx ^ dy\ 3B / dx dy\l 

+ L4(wa— wp)V^\ ^9" ^9"/ 2(mp + »2<T) Y da ^ da) j 

+ -7 ^ — ^ I 60:^:5 60;^ :5^ —3m 5-) — Ti ; r 

' |_ {ma — np)V 1 \ da da ^ dw/ 4(TOp+n(7) 

On putting a)=0, and therefore x=y = ^, this becomes 

(-P-"V=1)(45).,.=0. (19) 

Since A>\* for each of the equilibrium points for 0<(U<i, it can be 
readily seen that two of the roots of equation (5) are real and equal numerically 
but opposite in sign, and that the other two are conjugate imaginaries. Further, 
none of the roots of (5) is zero. It follows that none of the relations, 

jop=0; pp = — p; pp=±(tV— 1, 
can hold for p>l (2? a positive integer). Hence, from (19) we see that 

(1*^) =0. (20) 

Similarly, by differentiating the second, third, and fourth equations of (15) 
with regard to o, we obtain, respectively, 

* Moiilton, " Periodic Orbits," p. 159. 
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It follows then, that 

/a^\ ^0. /^\ ^0. /^\ ^ arbitrary =.. 
\ Old I M=o \ oa / 0=0 \Oa / a=o 

Hence, 

(1^) =c; and(l^) =-MC. (22) 

\Ca/ ai=0 \d(i)/ a=0 

If we differentiate equations (18) and the three corresponding equations in 
%, ti3, and Ui obtained from (15), and put o — O, we obtain in succession, by the 
aid of (20) and (21) 

/aX\ 1 r 3mim^—2)B 3mB 1^ 

\9«V«='o~ 2p + a'\/^L2(w(T— wp) V^ mp+naj. 
/d\\ _ 1 r 3mim^—2)B _ 3mB 1 ^ 

\8oV»=o~ 2p— (7V^^L2(»wr— «p) V^ mp+McrJ ' 

/aX\ _ J^ r 3n(w"— 2)i? 3mB 1^ 

\3a)*/ffl=o 3p L 2{ma — np) mp+na J ' 

/d%\ __ J^ r 3n(m''— 2)g 3mB 1 , 

\96)*/ffl=o p L 2 (ma — np) mp + naj 

On reduction, then, we obtain 

/d^x\ ^ ^ r 3w(w^— 2).gg _ n{m^—2)B 
\3cjV"=o L(w(r — t«p) (4p*+(T^) (m<T— wp)p 



(23) 



12mBp imB 1 

(wp + »2(T) (4p^ + <T*) (wp + n(T)pJ' 



/ a^\ ^ _ ^ X 6?^(w'— 2)^p __ 2n{m^—2)B 
\3o*/<.>=o L(m<T — wp) (4p^+cr*) (ma — «p)p 



6nBa 2mB 1 



(24) 



(wp+«(T) (4p^-|-(j*) (wp + v«cr)p 
By repeating the process and applying (20) and (21) at each step, we can 

/9*'M-\ 

find in succession the values of I ->r— i I (* = 1, . . . . , 4 ; k = 3, . . . . oo ) , and 

thence we can readily find I ,;^-ir and (,r-^ . We notice from (22) and 

\9w*/«=o Xoa^/a^o 

(24) that the first and second partial derivatives carry the arbitrary parameter 
c to the first and second powers respectively. We see further that, at each 
step in the differentiating, the right-hand members of the equations are homo- 
geneous in the orders of the partial derivatives in each term. It follows then 
29 
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that ( ^=-^1 and, therefore, also ( ^r-^ ) andf^^;-^) each carry a factor c* 

after the values of the partial derivatives of Orders lower than k have been 
substituted in the right-hand members. Henc6, the terms of the expansion in 
(16) carry c arid o as factors to the same power. On substituting the values 
of the partial derivatives in Maclaurin's expansions for x and y, and replacing 
6) by its value e"*", we obtain a set of solution^ of equations (8) in the form 



_„, , 1 r 3m(m^— 2)B(T n(m''-^2)i 12mBp 

P*' -X~ 1 -— — ( 1 ^ — • — ■ = ■ 

2lL{ma — np) {4:p^+a^) (wcr — np)p (mp+na) (ip^+a^) 



x= ce~'" + 

{mp+na)pj 3! 
__ _„, mV 6OT(m'— 2)gp 2n{m^~2)B 6nBa 



2!L(w(T — wp) (4p^+(T^) (wcr— kp)p (mp+*2a) (4p^+ff^) 
2mB 
{mp-\-na)p 



+ -7-^T^^^-l cH-^<"- ^ [ . . . . ] (?e-^<"+ 



(25) 



(B) Solutions in Powers of e"'"'". 

If, in equations (8), we were to transforni our independent variable by 
writiag a = 6+''*, we could build up a second set of solutions of form similar to 
(25), by a process exactly parallel to that used in section (A). We shall show, 
however, that this new set of solutions can be obtained directly from solutions 
(25) by changing the sign of t throughout, and changing the sign of y in the 
result. 

The first two equations of (1), with z dropped from the right-hand mem- 
bers, can be written in the form 

x"-2y'=F,{x,y''), 

y" + 2x' = yF,{x,y'). 
If we suppose the initial conditions are 

a=(0)=a, x'{0)=a„, y{0)=^, y'{0)=^^, (27) 

then the solutions of (26) have the form 

x=f,{t), x'=^At), y^W), y'=^2{t)' (28) 

If in equations (26) we put x=x, y — — yi, and t = — r, we get equations 
of identically the same form in x, r„ and t as (26) are in x, y, and t. These 
equations are 

x—2n=Fi{x,Yt), 1 

n + 2x^^F,{x,yi'), r 



. (26) 
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where the dot denotes differentiation with regard to r. If we impose the same 
initial conditions in these new variables, viz. : 

a;(0)=a, x{0)=a„ >7(0)=/?, v{0)=^^, (30) 

then the solutions of (29) are 

!V=fi{t), a;=^(T), v=fA'P), '7=^2(1'); (-31) 

that is X, X, ri, and yi are the same functions of r as x, x', y and y' were of t 
before. Bnt the initial conditions (30) are the same as 

a;(0)=a, x'{Q)^~a„ yo=-^, y'{0)=(3^. (32) 

Again, the solutions (31) are the same as 

x=f,i-t), x' = -^,{-t), y=-f,{-t), y'=^,{-t). 

It is readily seen, therefore, that for initial conditions (32) equations (8) 
admit of a solution which can be derived from (25) simply by changing the 
sign of t throughout and changing the sign of y in the result. These solutions, 
therefore, have the form 

_ ^j 1 r 3m(TO'— 2)gg n{m'—2)B 12wgp 

~ ^ 2\L{ma — np) (4p^+(T^) {ma — ■wp)p (mp+wcr) (4p^+o^) 

(mp + wcr)pj ^ 3! 

.^^0 6"*+^ r 6n(m^-2)^p _ 2n{m^-2)B 

^ 2! L('/n(T— np) (4p^+o^) (mcr— wp)p (mp+wcr)(4p^+(T^) 

(wp+n<T)p J ^ 3! 

V. Proof of the Convergence of the Solutions. 

It is necessary to prove the convergence of the series (25), which we have 
deduced from equations (8), or their transformed equivalents (9). This we 
shall do by a method analogous to that used by Picard.* Write equations (9) 
in the following form 



4m5 1 1 , ]^3g3p*^ 

(mp + na)pj 3! 



6nBa 



(33) 



9^-^ -cTV-lwi=.fff>+Hf + . . . . 



du, 



-p"^ +aV-l%=fff +flf + . . . . 

C/6) 

■P"|^+PM4 =i7f+Hf>+.... 



[ (34) 



♦Picard, "Traite D' Analyse," Vol. Ill, Chap. I, § 12, 
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The coefficients of any of the partial derivatives of Ui ii = l, . . . ., 4) for 6)=0 
are respectively of the form 

2?(-p)-ffV=l, p(-p)-(-(tV=1), pi-p)-p, p{~p)-{-p), (35) 

all of which are, in absolute value, greater than s{p — 1) where s is a real 
quantity. In this case e is a fixed quantity smaller than the absolute value of 
the smallest of ±p and ±<tV — 1. 

It has been shown by Moulton* that the expansions in the right-hand mem- 
bers of (14) converge within and on the boundary of a circle of radius a (a>0) 
about each of the equilibrium points ; the value of a depending upon which of 
the points (a), (&), or (c) is being considered. 

Let M be the maximum modulus of the expressions in the right members 
of (34) for all values of the variables %, u^, Ms, Ui, in this circle of radius a. 

Consider the comparison set of differential equations 



{"^-^1 



^^ _M-M^^l±^^^+^3, (. = 1, . . . .,4). (36) 



Vi + V2+Vs + Vi 



It is evident that equations (36) dominate (34) ; and it can readily be shown 
that the solutions of (36) dominate the solutions of (34). The terms on the 
right of (36) are all positive, and it is obvious from the method by which the 
solutions are built up that all the terms of the solution of (36) are positive. 
In building up the solutions of (34) and (36) as power series in a, whose 
coefficients are the values of the successive derivatives of the w,- and the Vi, re. 
spectively, for u=0, we see by (35) that the absolute value of the coefficient of 
any partial derivative on the left obtained from (34) is greater than the 
coefficient of the corresponding partial derivative on the left obtained from 
(36). But each term on the right-hand side of (36) dominates the correspond- 
ing term on the right of (34), and therefore any partial derivative of the right 
members of (36) is greater than the absolute value of the corresponding par- 
tial derivative of the right members of (34) . It follows, therefore, that the 
values of the successive partial derivatives obtained from (36) are greater 
respectively than the absolute values of the corresponding partial derivatives 
obtained from (34). Hence, each term in the solutions of (36) is greater than 
the absolute value of the corresponding term in the solutions of (34) ; or the 
solutions of (36) dominate those of (34). 

* Moulton, "Periodic Orbits," p. 154. 
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It remains to be shown that the solutions of (36) are convergent for a 
suifieiently small. From the symmetry of (36) in the v), {i = l, . . . ., 4), all 
the Vi are equak We can, therefore, replace equations (36) by a single equa- 
tion in one variable, viz. : 



{"i^-^] = 



M — M — , 



a 



which reduces to the form 



U;5 v= T— , where M'= : (37) 

Cut a—Av as 



whence, 

3w dv M'dv 



6) V a+{M'—4:)v 
On integrating this equation, we have 



log at = log jp- 

[a+{M'—4:)v]«^* 

or c(o = 7 r—-^ , where x = 



h^^y 



M'— 4 



Therefore, 



V — CO) ( a" + xa""^ — v+ .... 1=0. 



X 



From this it follows, by the theory of implicit functions, that v can be expressed 
as a power series in a, vanishing with a, and convergent for o sufficiently small. 
Thus, we see that the solutions of (36) converge for u sufficiently small; and it 
has been shown that the solutions of (36) dominate the solutions of (34). 
Hence the solutions of (34) are convergent for a sufficiently small, that is for 
t sufficiently large. 

VI. Alternative Method of Constructing the Solutions. 

The solutions of equations (8) may be built up by a process quite different 
from that used in § IV. 

If, in equations (2), we make the transformation x=X6, y=yB, the first 
two equations, considered for s—Q, give, on dividing through by e, 

x"-2y'={l + 2A)io+iBi-2a?^y-']s+2G{2a?-Zxy-'^^-^...., | 
y"+2x'=(l-A)y+ZBxye+iCy[~'^x^+y^]s''+ J 
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Since, as was pointed out in § V, the right-hand members of the original equa- 
tions converge within and on the boundary of a circle of radius a (a>0) about 
each of the equilibrium points, it follows that the right-hand members of (38) 
converge for e sufficiently small (0<e<l). Such a system of differential 
equations, with constant coefficients on the left, can be solved for x and y as 
power series in s, of the form 



00= I,xAt)e^, 



y= '^Viii)^, 

j=0 



(39) 



where the coefficients Xj and y^ are to be determined. Since (39) are the solu- 
tions of (38), the latter must be satisfied identically when x and y are replaced 
by the power series in (38). When we make this substitution, from terms 
independent of e, we have 

x'o'-2y'o-{l + 2A)Xo=0, 1 

y',' + 2x',-{l- A)yo=0. J 

The general solution of these equations is 

Xo=K,e''^=''*+K,e-''^=^*+K,e'"+K,e-''\ ] 

where the Ki and Li are constants of integration, the Ki (i=l, . . . ., 4) being 
arbitrary, that is dependent upon the initial conditions imposed which are 
arbitrary, and the L^ depending on the K^, the relations between them being* 

(1^+1 + 2^ 



p^-l-2A 
^'- 2p 



K,^ 


:V^ 


-lnK,= 


K^ 
-K,^'' 


Ks = 




mKs = 


K, 
K.^^- 



(42) 



Since the Ki are arbitrary, and since we are seeking solutions which vanish 
as t approaches infinity, we choose Ki = K2= Kg^O, and, therefore, also 
L^=L2=L3=0. Equations (40) then have a particular solution of the form 

Xo^K.e-"', y,=L^e-<'\ (43) 

Let us impose the initial conditions that x~b at ^=0. Since 

'K{t)=X,{t) +x^{t)s + x^{t)^+ ...., 

then 

Xo{0)=b, x,iOy=0, (i-l, ....<x>). (44) 

* Moulton, " Periodic Orbits," p. 159, Eqns. (28) . 
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On imposing these initial conditions and using (42), we see that (43) takes 
the form 

Xo^be-"*, ya^—mbe-"*. (45) 

In finding successively the values of Xj and yj (;/ = l, . . . .co ) it is more con- 
venient to use the normal form of equations (26). 
As we saw in § 2, the transformations 



X =«i + M2 + M3 + 'M4, 

a^' = 0'V— l(«i— M2) +p(W3— W4), 

y' = — na{Ui-\-Ui) +mp{us + u^), 
change equations (26) into the normal form 

m [.....] 6 \ . . . .\s 



(46) 



u'l — (tV — 1mi= 
«2+crV— 1%= — 



2 {ma — Mp)V — 1 2(wp+w<T) 
m[. . , .]e \ . . . .\s 



2{ma — «p)v — 1 2{mp-\-na) ' 
—n[ ]e , j \e 



^' P^' ~ 2{ma—np) '^ 2{mp+naj ' 
, n[ ]e \ |e 

Ui + pUi = ^r-r r \- 



2{ma — np) 2{mp + n(y) ' 



(47) 



where 



(48) 



[ ] = |Z?[— 2a;2+«/2] +2C[2.x^—3xy^]e+ terms in e^, e^ , 

\ \=3B\xy\ + iGy\—4:X^ + y^\s+ terms in e', e^ 

These equations have solutions of the form 

M,= iwi'V, (*=1, ....,4). 

1=0 

To determine the «p^ we substitute (49) in (47) and equate coefficients of like 
powers of e on both sides of the resulting equations. From the terms inde- 
pendent of £, by applying (46) and (44) we get the value of Xo and «/o as given 
in (45). Since, by (44), 

a;o(0) =&, and Xj{0) =0, (i = l, 00 ), and Xj= 2 mP, (49) 

it follows that 



<>(0)+«r(o)+Mf (0)+t*f>(0)=6, I 

u^HO)+uiHO) + niHO)+u^HO)=0, i=(l, ....a,). J 
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On equating coefficients of the first power of e, we have 



3mB[-24+yg] ^Bx,y^ 



A{ma — «p) V — 1 2{mp+na) ' 
4(m<r — np)v — 1 2{mp+na) 






3nB[-2xl+yl] ZBx,y, 

4(w<T— «p) 2{mp-\-n<j) ' 

3nB[-2xl+yl-\ 3Bx,y, 

4:{ma—np) 2(mp + n<T) ' 



(51) 



After substituting for Xq and t/o their values from (33) we have 

[ 4(m<T— wp)V^ 2(mp + n(T)J '' ' 

[ 4(m<T— «p)V^ 2(mp+«(T)J '' 

= ( ^-B(m^-2) _ SmB\^,^_, 

I 4(m(T— wp) 2{mp + na) J 



mJ^^'-P«^" 



«r^'+pM?^ 



(mff — wp) 2(mp + w(T) 

-f 3w^(w''— 2) 3w^ 

I 4(mo' — wp) 2(mp + «<T) 



(52) 



where the Mg^ («=1, ,4) denote the expressions in the brackets. The 

superscript and the first subscript on the M^^ are the same as those on the 
corresponding uf^, while the second subscript denotes that they are coefficients 
of e-^"*. On integrating (52), we get 



<)=:cf)e'V-u -^il'^' 



2p + 



(TV — 1 



MW=c«e-'^-"- 



4" = 0^1)6"' 



itO) = c«e-'" 



M(p&^ 



2p 



-~?-^=e-^'"=c("e-''^-^'— cfg^fe^e-^"', 



M">6^ 



3p 



- e-^'"=C8«>e'"— dg^fo^e-^"', 



g-2p*_^(l)g-pt_^(l)^2g-2„« 



(53) 



If then we put cJ^>=C2")=cp>=0 and choose cf^ so that the initial conditions 
(50) are satisfied, we find that 
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ci'^=b'[di^+di^+4v+dgn. 

Since 

we see that 



(54) 



2 
2 



ooi=[(d\i'>+dg>+di^+dg->)e-'"^{d\l^^d<^ + dg^+di^)e-'^']b 

= [Z>0)e-p«_j9a)e-2p'] b 

y,= [-mid[l>+d^l> + dil> + dl^)e-''*-\nV^=i{d\^-d'^2 

+ midi^—dil^) |e-''^«]6^^ [— wD^i^e-"'— Z>|'>e-^'"]&^ 



(55) 



n'2 



f +<tV^:3w(2)^_^ 



From the coefficients of e^ above, we see that the differential equations 
which define the «f^ (t=l, ....,4) are 

(2)' ■\/'~T (2)_ 3mB{—2XoXi+y^yj^) 3B {Xoyi + XiPo) 

2(m(T— wp)y — 1 2{mp-\-na) 

^ mC{24-3x^ ^ 3C{4=4yo-yl) 
{ma — wp) V — 1 4(wp + «(T) 

3mB {—2xQX^+y^yT) 3B {x^y^^ x^y^) 
2{m<s — wp) V— i 2{mp+na) 

mCi24-Sxoy^ ^ 3G(4:xlyo-yl) 

(ma — wp) V — 1 4:{mp+na) I /qq\ 

^(2)-_ (2) ^ _ 3nB{—2xoyo+yoyi) 3B (Xpy^-^ x,y^) 

^ y^ 2(m0—np) 2{mp+na) 

nGi2x^,-3x,yl) 3C{4xlyo-yl) 
ma — np 4(wp+W(T) ' 

(2)' , ,,(2) _ , SnB{—2xoyo+yoyi) 3B {x^y^^ x^y^) 

u, -tpu, --t 2{ma—np) '^ 2{mp+na) 

nG{2xl^3xoyl) SCjAxlyo-yl) 
ma — np 4:{mp+na) 

On substituting for Xq, y^, x^, y^ their values from (45) and (55) these equa- 
tions take the form 



ufy—a^/—luf^^=Mf^b^e-^<"+M^^^b^e-^''\ 
uf^'—pup =Mf^Ve-^<"+Mf^h^e-^''\ 



(57) 



80 
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where 



2'{ma — wp) V — 1 2{mp+na) 

3m.g(w^— 2)D|^> 6mBDi^^ ^ 

2 {ma — wp)V — 1 2{mp-\-na) 

ZnB{m^—2)I)f^ 6mBD^^^ . 

2(m<T — ■wp) 2{mf)+na) ' 

3»g(w^— 2)Dq> 6mBDi'^ ^ 

2{mcf — np) 2(wp+na) ' 

3mB{2Dl^^ + mM^) ZB {mBf^ —W^^) mC{3m^—2) 






2{m(j — Mp)V — 1 2{mp-\-na) (ma — np) V — 1 

3wC(w'— 4) _ 
4(TOp+»«<y) ' 



+ 



_ 3wg(2D|^)+wZ)|'>) 3.g(wDp— J^^)) wC(3m^— 2) 

ilx 23 — „ , ^ _ / — ^ , V ~l~ 



+ 



2{ma — wp)V — 1 2{mp-\-na) {ma — np)'\/ — 1 

3m(7(w^— 4) 
4(wp+wo') 

_ 3nB{2D^^^+mW^^) 3.g(wZ><'>— ■Pi^>) nG{3m'—2) 
^ ~ 2{ma — np) 2{mp-{-na) ma — np 

3mC(m^— 4) 

4(wp + «(T) 

jl£(2^^ 3wg(2Z)Q> + m5ffl) 3ig(wI>(^>-.D|^)) nC(3w^-2) 
*^ 2(w(T — np) 2{mp+na) ma — np 

3mC(w'— 4) 



4:{mp+na) 



The solutions of (57) have the form 

<) = c\^'e''''-'*—d^^b^e-^'"—dlf¥e-^''\ 
Mf ) = 4^^e-'''^-'*—d^fb^e~^''*—dg^Ve-^'", 

wf = cth-''*—difbh-'''—d^i^b'e-"'*, 



where 




2p+av—l 2p— (TV— 1 3p 
3pH-oV— 1 3p— (TV— 1 4p 


^(2)_ ^ii> 

^"•- 2p 



(58) 



(59) 



h (60) 
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If we put cp> = 4^> = c^^^=0 and determine cf^ to satisfy the initial conditions 
(50), we have 

On substituting these values for cf^ (*=!, . . . ., 4) in (59), and writing 

2 4'^=Df) ; nV^id^f-d^f) +m{dif-dg^) =Bf\ 

2 4')=Df>; wV-l(c?g>-4r) +m{dg^-d'i^) =.Df\ 
w.e have, by (54), 

y^=[—m{Df^-\-Df'^)e-'"—Df>e-^''*—D'i'>e-''''}h\ J 

If we proceed in this way, we can build up in succession the values of Xj 
and t/,- for i=3, 4, . . . . oo . By induction we can get the form of the general 
term. From (45), (55), (62) we notice that the x^ and y^ are sums of powers 
of e"'", the highest power of e"'" occurring being j-\-l. The equations defining 
ul'^ (i=l, ,...,4) have in their right-hand members only sums of powers of 
e""', the lowest power being e"^"* and the highest e^^'+^^^K When we integrate 
these equations we will have terms in powers of e"*"* the lowest power being 
the first and the highest the (r+l)-th. These solutions will have the form 

ttl") =cf >e^'- [ 2 4''>e-^'"] ft'+S (i=l, . . . . , 4 ; 

;^-(tV^, — <tV^, p, — p). (63). 

If we put c^'') = c|''> = c|''^=:0, then to satisfy the initial conditions we must put 

4 v+1 

c('') = &''+i2 2 4''\ If then we put 

Dr=24^ (i=2,....,».+i) 1 

i=i I V (64) 

Dj''^=nV'^{dip-dif)+midip-d^^), (i = 2, . . . ., v + 1), . 
it follows that 

/=2 fc=2 

«'+l I'+l _ 

y^= [—m 2 Z)}''>e-'"— 2 Di''^e~'""]b''+\ 

;=2 i;=2 

When we substitute the results of (45), (55), (62), (65) in (39), we have the 
values of x and y which are solutions of (38). The x and y belonging to the 
physical problem as defined by (8) are s times the x and y, respectively, which 
we have just obtained from (38). On multiplying the values of each Xj and y^ 
by s and substituting in (39) , we see that the resulting expressions carry b and e 



(65) 
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as factors to the same power in each term. Hence be is equivalent to a single 
parameter and may be replaced by /?. Our solutions then may be written 

fc=2 h=i 

v+1 C+1 

+ [ 2 B'i'^e-"*— 2 Z>i''>e-'"'*]/3''+i+ . . . ., 

fc=2 fc=2 

*=2 J;=2 

I>+1 v+l_ 

— .... — [w 2 Di''>e-'"+ 2 Z)i''>e-'""]/3''+^— .... 

fc=2 fc=2 

If, instead of taking a;(0) =fe as our initial conditions, we had taken 

x{0) =b^D\'^bh—Df>V^—Df^¥i—.. . ., 

we would have had 

Xo(0)=6, x,{Q)=-D<fl,y+\ 

Then at the successive steps the constants of integration 

cf) = cf = c|=....c«>=....-0. 

The u^'^ (i=l, . . . ., 4| ^=0, . . . .oo ), instead of being expressed as sums of 
powers in e~^\ would each be expressed as a single term in 6""^^+^)"*. The solu- 
tions (66) would then have the form 

x= e-*'i3-j5|i>e-^^'^2— Dp >e-»'"/3«- .... -Z)(';\e-<''+i>'"/3''+i- . . . . , 1 

y = _ we-'"/?-5«)e-^'/?^-5<^>e-«'"/3^- ^^'^le-c+^w'/J'+i- J ^ 

If We compute the values of the coefficients in these expansions in (68) and write 
G for /8, we find that solutions (68) are identical with solutions (25) of § IV. 

By a method exactly analogous to that just used we could build up a set 
of solutions arranged in ascending powers of e+'", and with a proper choice of 
initial conditions it could be shown that they were identical with solutions (33) 
of § IV. 

Vll. Properties of th6 Orbits. 

When we consider the solutions of the differential equations as given by 
equations (25) of § IV, which are in the form of power series in e""*, we see 
that the values of x and y continually decrease as t becomes greater and 
greater, and finally x and y approach the value zero as t becomes infinitely 
great. Since x and y are the coordinates of the infinitesimal body referred to 
an equilibrium point as origin, it follows that it will approach nearer and 
nearer to one of the equilibrium points as t increases. Such orbits are said to 
be asymptotic to these points. 

Similarly, we see that if the infinitesimal body were moving on one of the 
orbits which are given by equations (33) of § IV, where x and y are expressed 
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in the form of power series in e"*""*, it would approach one of the equilibrium 
points when t became infinitely large and negative. In other words, if we 
imagine the infinitesimal body to be placed at one of the equilibrium points, it 
Would gradually leave it on one of these orbits, requiring, however, an infinite 
time to describe the first small part of the orbit. 

(A) Meaning of the Parameter c. 
When the masses of the finite bodies are given, the only arbitrary in solu- 

^-) . If c is fixed, the x and «/ are determined 

0(0/ «=0 

uniquely by (25) for all values of t so large that the convergence of the solu- 
tions as power series holds. Now the slope of the orbit at any time t, for t 
sufficiently large, is given by 

dy_ftl _ wpce-^'+2pJ!|^>c^e-^'"+3pZ)f c^e-^^'-f • • • • /go-, 

dx~ 1% ~ — pce-'" + 2pi)«)c2e-^''*+3pi)f>c«e-^'"+ . . . . ' ^ ' 

= —m-\- a power series in e""' beginning with the term in e~''\ There- 
fore, 

(^) =-m. (70) 

\dx/ 1=«, 

Hence the direction of approach to the equilibrium point is independent of the 
value of the arbitrary parameter c. 

The position of the body in its orbit at the initial time, ^—0, is given as a 
power series in c with constant coefficients ; so that if c is fixed the position of 
the body in its orbit at the initial time is determined. It follows therefore 
that the direction of approach to the equilibrium points is independent of the 
initial position of the infinitesimal body in its orbit at the initial time. Since 
c is arbitrary, nothing of generality will be lost in the actual construction of 
such orbits if we take c equal to unity, and the numerical computation will be 
simplified. This we do in the numerical computation of an orbit in § VIII. 

(B) The Number and Position of the Asymptotic Orbits. 
In equations (42) of § VI the quantity m is defined 

w= c . 

2p 

It will be shown later* that for each of the equilibrium points (a), (&), (c) the 
quantity m is negative for all values of /u, where 0<|«<i. Therefore, in the 

ease of the oirbits given by (25), it follows that (--^) = a positive quantity. 

\dx/ t=K 

Hence equations (25) represent two orbits, one in the first and one in the third 

♦See §VII (C). 
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quadrant in the neighborhood of each of the equilibrium points. For any pre- 
assigned valuer of ft these two orbits are equally inclined to the a;-axis in the 
neighborhood of the origin. 

Similarly, it may be shown that equations (33) represent two orbits 
leaving each of the equilibrium points, one in the second and one in the fourth 
quadrant. For any given value of f^, these orbits are, in the neighborhood of 
the origin, the images by reflection on the a;-axis of the orbits given by (25). 

Near the points of equilibrium, that is for t very large, the first terms in 
the expansion of (25) are the most important and determine the sign of the 
right-hand members. It follows therefore from (25), since m is negative, that 
if c is positive x and y are both positive, but if c is negative x and y are both 
negative. Hence for positive values of c the orbit is in the first quadrant, and 
for negative values of c the orbit is in the third quadrant in the neighborhood 
of the equilibrium points. 

Similarly, when t is very large and negative, the first terms in the expan- 
sion in (20) determine the sign of the right-hand members. For c positive 
X is positive and y negative, and for c negative « is negative and y positive. 
Hence for positive values of c the orbit is in the fourth quadrant, and for 
negative values of c the orbit is in the first quadrant. 

The value of -^ for a point near the origin is given by (69). If we take 

the second derivative 




we find that 



d^y 
dx^ 



■■ — 2[Z>^^)+toD|^^] + a power series in e'"*. 



(72 \ 
— % ) is independent of c, and therefore does not change 

sign with c. Therefore, near the equilibrium points, the orbits in the first and 
third quadrants lie on the same side of their common tangent line 2/ + ma; = 0. 
-Similarly, it can be shown that the orbits in the second and fourth quadrants 
lie upon the same side of their common tangent, y — ma/=0. 

(C) Variation in the Direction of Approach as (i Varies from Zero to 1/2, 

Now it has been shown* that if fj and r^ denote the distances of the 
infinitesimal body from the finite masses 1 — (i and ^ respectively, the values 
of rx and r^ for the equilibrium points can be expressed, for ^ suflficiently small, 

* Moulton, " Introduction to Celestial Mechanics " (New Edition) , § 158. 



Waeben: a Glass of Asymptotic Orbits in the Problem of Three Bodies. 243 



as convergent power series in (i^^^ in the case of the points (a) and (&), and in 
(I in the case of the point (c). 
For the point (a), 



--(f) 

ri=l-\-r. 



3 \3/ 9 V3/ 



For the point (6), 



\S/3 



For the point (c), 



_^ 7 23X7^ 3 



12* 



= ^2—1, 






M 



m: 



p^— 1— 2^ 
2p ■ 



(72) 



we have seen further that — m gives the value of the tangent made by the 
curve with the positive «-axis in the neighborhood of one of the equilibrium 
points. 

The following table has been constructed to show how the elements r^, r^, A, 
and — m vary as the ratio of the finite masses changes from zero to one. 



M 


M=0 


/x=0.5 


Point 


(a) 


(&) 


(0) 


(0) 


(6) 


(0) 


r2 


<r 


<r 


2 w" 


0.63273 


0.4308 


1.700 


ri 


=1 


=1 


-^^ 


1.63273 


0.5692 


0.700 


A 


=4 


=4 


1+lfi^l+e 


2.0886 


8.9654 


1.5595 


P 


= 2.5083 


= 2.5083 


V2e 


1.5564 


4.0305 


1.1469 


— m 


= 0.56224 


= 0.56224 


= 00 


0.88498 


0.3324 


1.2232 


4> 


= 29° 20' 


= 29° 20' 


= 90° 


41° 30' 


18° 23' 


50° 44' 



From the above table we see that as fj. increases from zero to 1/2 
for (a), ^ increases from 29° 20' to 41° 30' approximately; 
for (b), ^ decreases from 29° 20' to 18°23' approximately; 
for ((?), ^ decreases from 90° 00' to 50°44' approximately. 
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VIII. Continuation of the Orbits beyond the Range of Convergence of the 
Solutions of the Differential Equations. 

It has been pointed out that the expansions given in (25) and (33) define 
orbits asymptotic to the equilibrium points for t sufficiently large. It remains 
to be shown what becomes of the infinitesimal body for smaller values of t. 

(A) Method of Mechanical Quadratures. 
If the X- and ^/-coordinates of the moving body be computed for a sufficient 
number of equidistant values of the time then the coordinates of the body for 
the next equidistant value of the time may be computed by the method of 
mechanical quadratures. If, for example, we know the values of a/i«/i, ^^22/2, 
'"^sysj *4^4> %2/6 3.t five equidistant values of the time ti, t^, t^ft^, t^ then we shall 
show how the value of cce and t/g for the time te can be computed, where 

tg ti = t^ ^4= .... ^^2 *!• 

Form a table of the values of x and their successive differences for the 
successive values of the time. Form also similar tables for y, x', y', F^, F^ 
where from (1), 
dx 

dt -"^ 



d^x dx' _ , ,\ o , , (1 — u){x — Xi) u{x—X2) 

^=^y' 
dt ^ 

d'y _ dy' _p^^^^ ^,^ y^=_2a,>+y^ Ik:^ _ M 



(73) 



dt^ ~ dt -*^-^-^^' -- ' ^ ^3 ,.| • 

It is necessary to have computed the values of these quantities at sufficient 
dates that some order of differences obtainable from them will be small. Sup- 
pose that such tables of values have been set up for each of the quantities 
^, y, *') y'j ^2» t^i fo^" ^^^ dates small ti, t^, tg, . . . ., t^, and suppose that in these 
tables the fourth differences are small and approximately constant. In the F^ 
and Fi tables we assume the next fourth difference about equal to the previous 
ones, and from the tables compute F^"+^) and Fl'^'^'^K Then the values of the 
next first differences of x' and y' can be readily computed by means of the 
formulae : 

Ar»a;'=(f„+i-^„)[Fr+»-iAiF^"+^>-^A2F^"+^)-^A3F<»+^> 

36^*^^ J' 



-3^A.F<->....]._ 



(74) 
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Having found A[''+^'>x' and Ai''+^'>y' we can at once compute x'{t„+i) and y'{t„+^), 
and complete the table of differences for this date. We can then compute 
^(n+i)^ and A{"+^)t/ from x'it^+i), y'it^+i) and their successive differences, by 
formulae exactly similar in form to (74), namely, 



36 ' 



(75) 



■] 

Ar'^y= (^„+i-^j[/<»+^>-| Aii/'<"+^)- ^ A^'(»+'>- ^ A3/<"+^) 

L A ^'(»+l). ,,'.1. 

36 *^ J'. 

whence we obtain at once x{tn+i) and y{tn+i). 

Having found the values of x, y, x', y' at the date t^+i we compute F^''+') 
and i?'<"+^) by the second and fourth relations of (73). If the results so ob- 
tained agree with the results we already have for F|"+^) and F^"+^>, we assume 
that the fourth differences which we guessed are correct. But if the computed 
values of F^^+i) and ^^"+1) differ from the assumed values we replace the latter 
by the computed values, and repeat the process as before from this point on. 
We keep on repeating this process, which generally has to be done only once, 
until the computed value is the same as the value from which it is computed. 
Having thus found the coordinates of the body for the date i„+x, we assume 
another set of fourth differences for F^ and F4 for the date t^+i and repeat the 
process. Thus any number of points on the orbit may be determined', and the 
orbit can thus be continued beyond the range of convergence of the solutions 
of the differential equations. 

(B) Jacobi's Constant. 

, The question arises how we know that the orbit so continued is really the 
orbit on which the infinitesimal body would move when forming a part of the 
system in question. 

It has been shown * that in the case of the infinitesimal body 

r=^2+ 2_^2(iz::^ + ^_c, (76) 

or 

C =«^ + /+ ^i^— ^ + -^ -x''-y'\ {11) 

" rx r^ 

This gives us a relation that must hold between the coordinates and the x- and 

^/-components of the velocity of the infinitesimal body at all positions of its 

Orbit. 

*Moulton, "Introduction to Celestial Mechanics" (New Edition), § 153. 

31 
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The constant G, which is called Jacobi's Constant, can be computed from 
the initial configuration of the system. In the present case the value of C 
would be the value of the right-hand member of (77) when the body is at rest 
at the equilibrium point in question. At such points x'=y'=x='y=0; r^, and r^ 
are readily determined from (72). 

Formula (77) may then be used as a check on the computation at each 
step, that is, for each new date. If the integral will not check, when the com- 
putation has been accurately made, it is usually necessary to divide the time- 
interval, and take shorter steps. It then becomes necessary to interpolate 
values of the coordinates for intermediate values of the time. After obtaining 
the intermediate values for F^ and F^ by Lagrange's Interpolation Formula, 
the corresponding OD, x', y, y' can be computed out of these by the process 
described. Thus the orbit of the infinitesimal body can be continued as far as 
may be desired. 

(C) Specific Form of the Solutions for |«=0.02. 

"We proceed now to assign to ^ the arbitrary value ;m=0.02, and find the 
specific form of the solutions and the shape of the orbit. We will consider 
the orbits in the neighborhood of the equilibrium point {a), and consider those 
given by equations (25), where the expansion is made in powers of e""'. 

By (72) we have r^ for the point (o) given by 



' \3/ ^3\3/ 9 



On neglecting all after the third term in the expansion, we find, for /a =0.0 2 
that r2=:0.1993. But because of the neglected terms in (72) this value of r^ 
does not satisfy with sufficient accuracy the quintic equation of which (72) is 
the solution, namely, 

rl+ (3-w)r|-|- {Z'-2n)f^—^t\-2iir^-ii=0. (78) 

If we start with the value 1*2=0 -1993 and apply the method of differential cor- 
rections we find f2=0- 200078, a value which makes the left member of (78) 
differ from zero by approximately 0.00000005. Hence, for the point (a) we 

take 

r2=0. 200078, rj-l. 200078; (79) 

whence, by (2), (5), and (42) we find in succession 

A= 3.064095; p=2. 098701; <t= 1.827794; 
m = — 0.648885; w=2. 863898; ^ = 12.952988; 

mp + na= 3.872514; ma—np= 7.196426; 

4pH<T' =20.958701; m^ —2 =—1.578948. 

On substituting these values in equations (25) and putting c=+l,* we have 

* See § VII (A) . 
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for the equation of the orbit in the first quadrant in the neighborhood of the 

origin, 

^^ g-2.09870H_2 , 944651 e;-*-i»™2t ^ | 

2/ = 0. 648885 e-'-o'^^'+0. 0251562 e-*-'*^""'' J 

For the infinitesimal body at rest at the equilibrium point (a), Jacobi's constant 
C has the value 



6. = ,;^ + ^^+2(l-^ ^ 



:3. 225734. 



Hence, at all points on the orbits asymptotic to (a), Jacobi's constant must 
have the value 3 . 225734. 

The values of x, y, «', and y', were computed for the times t=Z, 2.75, 
2.50, 2.25, and 2.00 and the tables of values of x, y, x', y', 1\, F^ and their 
differences built up. On proceeding as outlined in subsection (A) of this article, 
the table of values (on p. 248) was obtained. For the first few computations the 
interval of time taken was one-fourth of the unit of time, and afterwards one- 
eighth ; but in the table are given only the values corresponding to intervals 
of one-fourth, except in the case of a few dates near which one of the elements 
changed sign. 

A drawing of the orbit represented by the data in the following table is 
appended. 




t 


X 


y 


X' 


y' 


3.00 


1.18191 


0.00120 


—0.00383 


—0.00251 


2.75 


1.18316 


0.00202 


—0.00642 


—0.00424 


2.50 


1.18526 


0.00342 


—0.01071 


—0.00717 


2.25 


1.18874 


0.00577 


—0.01773 


—0.01210 


2.00 


1.19447 


0.00974 


—0.02896 


—0.02038 


1.75 


1.20373 


0.01643 


-0.04633 


—0.03434 


1.50 


1.21830 


0.02766 


-0.07175 


—0.05741 


1.25 


1.24034 


0.04631 


—0.10604 


—0.09469 


1.00 


1.27190 


0.07674 


—0.14716 


—0.15243 


0.75 


1.31392 


0.12474 


-0.18816 


-0.23661 


0.50 


1.36500 


0.19749 


—0.21733 


—0.35027 


0.25 


1.42037 


0.30217 


—0.21993 


—0.49123 


0.00 


1.47154 


0.44471 


—0.18149 


-0.65126 


-0.25 


1.50674 


0.62827 


—0.09068 


—0.81669 


—0.50 


1.51198 


0.85206 


0.05873 


—0.96995 


—0.75 


1.47247 


1.11060 


0.26679 


—1.09141 


-1.00 


1.3742 


1.3935 


0.5269 


—1.1614 


—1.25 


1,2057 


1.6856 


0.8262 


-1.1621 


—1.50 


0.9594 


1.9677 


1.1459 


—1.0796 


-1.75 


0.6331 


2.2176 


1.4628 


—0.9050 


—2.00 


0.2305 


2.4122 


1.7508 


—0.6362 


—2.125 


0.0037 


2.. 4815 


1.8756 


—0.4681 


—2.25 


—0.2377 


2.5284 


1.9830 


—0.2780 


-2.375 


-0.4912 


2.5504 


2^0699 


-0.0705 


-2.50 


—0.7542 


2.5454 


2.1333 


0.1542 


—2.75 


—1 . 2957 


2.4468 


2.1794 


0.6420 


-3.00 


—1.8338 


2.2220 


2.1041 


1.1580 


—3.25 


-2.3369 


1.8680 


1.8979 


1.6706 


—3.50 


—2.7717 


1.3899 


1.5593 


2.1452 


-3.75 


—3.1060 


0.8015 


1.0961 


2.5468 


-4.00 


—3.3107 


0.1254 


0:5255 


2.8422 


—4.125 


-3,3567 


—0.2364 


0.2077 


2.9409 


—4.25 


-3.3619 


—0 . 6083 


—0.1265 


3.0027 


—4.50 


—3 . 2435 


—1.3629 


—0.8263 


3.0062 


—4.75 


—2.9482 


—2.0971 


—1.5344 


2.8394 


-5.00 


—2.4792 


-2.7680 


—2.2081 


2.4989 


-5.25 


-1.8507 


-3.3327 


-2.80.35 


1.9927 


-5.50 


-1.0875 


-3.7521 


-3.2789 


1.3399 


—5.75 


—0.2243 


-3.9929 


-3.5975 


0.5701 


—8.875 


0.2316 


—4.0383 


-3.6886 


0.1532 


—6.00 


0.6958 


—4.0307 


—3.7303 


—0.2778 


—6.25 


1.6237 


—3.8514 


—3.6580 


—1.1578 


-6.50 


2.5069 


—3.4534 


—3.3729 


—2.0189 


-6.75 


3.2927 


-2.8479 


—2.8803 


—3.9747 


-7.00 


3.9313 


-2.0593 


—2.1988 


—4.2632 


—7.25 


4.3790 


--1.1239 


—1.3597 




—7.50 


4.6016 


—0.0894 


—0.4070 




-7.625 


4.4446 


0.4475 


0.0947 
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